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ABSTRACT 

The equation of state and radiative opacities of partially ionized, strongly magnetized hydrogen plas- 
mas, presented in a previous paper [ApJ 585, 955 (2003)] for the magnetic field strengths 8 x 10^^ G 
< i? < 3 X lO^'^ G, are extended to the field strengths 3 x 10^^ G < i? < 10^^ G, relevant for magnetars. 
The first- and second-order thermodynamic functions and radiative opacities are calculated and tabu- 
lated for 5 X 10^ K < T < 4 X 10^ K in a wide range of densities. We show that bound-free transitions 
give an important contribution to the opacities in the considered range of B in the outer neutron-star 
atmosphere layers. Unlike the case of weaker fields, bound-bound transitions are unimportant. 
Subject headings: equation of state — magnetic fields — plasmas — stars: atmospheres — stars: neutron 



1. INTRODUCTION 

Neutron star atmospheres differ from the atmospheres 
of ordinary stars because of the high gravity and mag- 
netic fields. The majority of the kno wn pulsars possess 
magnetic fields B ~ lO^-lO^^-^ G l|Manchester et al.l 
l2003il . Most popular models of the soft gamma repeaters 
and anomalous X-ray pulsars assume they are magnetars^ 
neutron stars with superstrong magnetic fields B ~ 10^^- 
IQis G (e.g.,|Thompson 2003). 

The properties of matter are strongly modified if 
the electron cyclotron energy fUVcc = heB/rricC ex- 
ceeds 1 a.u. - i.e., the field strength B is higher than 
Bn = ml ce ^lh^ = 2.3505 x 10^ G (for a review, see 
lLaill200l() . iShibanov et"al] l|1992D presented the first 
detailed model of hydrogen atmospheres for neutron 
stars with strong magnetic fields. Lat er it was devel- 
oped and used by many authors (e.g., iHo fc Lai 120031: 
|Ozel"2003' 'Paee. Shibanov. & Zavlin"199i iPavlov et all 
^995; Zane et al. 2001: Zavlin ct al. 1995; and references 
therein). These works have played a valuable role in as- 
sessing the observed spectra of the neutron stars. Despite 
the recognition of the possibility that the absorption 
by atoms can significantly contribute to the opacities, 
the atoms were not included in the atmosphere model 
in a thermodynamically consistent manner in the listed 
works. 

Theoretical equation of state for partially ionized hy- 
drogen plasma with strong magnetic fields was developed 
by Potckliin. Chabricr. & Shibanov (1999). The nonper- 
turbative effects of atomic thermal motion across the 
magnetic field were taken into account and shown to 
be i mportant. Based on this th eory, in a previous pa- 
per ijPotekhin fc Chabried 1200.^ hereafter Paper I) we 
performed extensive calculations of the thermodynamic 

^ Also at the Isaac Newton Institute of Chile, St. Petersburg 
Branch, Russia 



functions and calculated monochromatic and Rosseland 
mean opacities, taking into account the bound-bound 
and bound- free transitions. These calculations were done 
at 11.9 < log i?[G] < 13.5. In addition, we reevaluated 
the free-free opacities, taking into account the motion of 
both interacting particles, electron and proton, in arbi- 
trary magnetic fields. 

Here we report an extension of the previous calcula- 
tions of the equation of state and opacities of a par- 
tially ionized hydrogen plasma to higher field strengths, 
13.5 < \ogB[G\ < 15. 

2. THEORETICAL MODEL 

A major complication when incorporating bound 
species in the models of strongly magnetized neutron- 
star atmospheres arises from the strong coupling between 
the center-of-mass motion of the atom and the internal 
atomic structure. If an atom rests without motion in a 
strong magnetic field, there are two distinct classes of its 
quantum states: at every value of the Landau quantum 
number n and the magnetic quantum number —s {n > 0, 
s > —n), there is one tightly bound state, with binding 
energy growing asymptotically (at B oo) as (InB)^, 
and an infinite scries of loosely bound states (numbered 
by V = 0, 1, . . .) with binding energi es approa ching the 
energ ies of a field-free H atom ("e.g.. iCanuto fc Ventural 
^23- The Landau number n is not a "good" quantum 
number, in the sense that the Coulomb interaction mixes 
different Landau orbitals of a free electron, but this num- 
bering is unambiguous and convenient at B ^ Bq. The 
atom is elongated: its size along the magnetic field B 
either decreases logarithmically (for the tightly bound 
states) or remains nearly constant (for the loosely bound 
states), while the transverse radius decreases as B~^/^. 
The radiative transition rates are different for the three 
basic polarizations: the linear polarization along the field 
and the two circular polarizations in the transverse plane. 
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This simplicity is destroyed when atomic motion is 
taken into account. The electric field, induced in the 
comoving frame of reference, breaks down the cylindri- 
cal symmetry. Separating the center-of-mass motion 
and ch oosing an appropriate gauge of the vect or po- 
tential JCor'kov fc Dzvaloshinskiilll967t [Potekhini ri994t 
iVincke. Le Dourneuf. &: Bavelll992(l . one comes to an ef- 
fective one-particle Schrodinger equation with an effec- 
tive Hamiltonian depending on the pseudomomentum K, 
which is a conserved quantity related to the center-of- 
mass motion in a magnetic field. If the component of K 
perpendicular to B, K±, is finite, then not only n, but 
also s is not a "good" quantum number; nevertheless, it 
is convenient to retain the quantum level numbering by 
n, s, and specified above. In t he adiaba tic approxima- 
tion, widely used in the past ( Gor'kov fc Dz valoshinskii 
119671 llpatova. Maslov. fc Suba shievI Il984|) . the trans- 
verse part of the wave function is postulated to be the 
same as for an electron without Coulomb interaction. We 
perform calculations without this approximation, using 
the technique described by Potckhin ( 199 4) for the dis- 
crete atomic spectrum and by Potckhin fc Pavlov! l)1997f) 
for the continuum. This technique uses the basis of Lan- 
dau functions, which represent the transverse parts of 
the wave functions of a free electron in a magnetic field 
and are the same for the solutions of the Schrodinger 
and Dirac equations. Therefore this technique can be 
used for the ordinary {huJcc < rricC^) as well as super- 
strong fields {hujcc > rricC^), provided the atomic binding 
energy <C rUcC^ , which is always the case &i B < IQ^^ G. 

An atom moving across the strong magnetic field ac- 
quires a constant dipole moment perpendicular to B 
and K. Those radiative transitions, which were dipole- 
forbidden for an atom at rest because of conservation 
of the z-projection of the angular momentum, become 
allowed and should be taken into account in the atmo- 
sphere models. The binding energies decrease with in- 
creasing K^. Asymptotically, at large Ki_, the binding 
energies tend to e^B/{cK±). 

The equation of state for partially ionized hydrogen 
in strong magnetic fields was constructed and discussed 
bv iPotekhin et al.l l)199i^ . The treatment is based on 
the free energy minimization. The free energy model is 
in essence a generali z ation of the nonmagnetic model of 
ISaumon fc Chabrieil lITMl IT991 to the case of a strong 
magnetic field. We consider a plasma composed of Np 
protons, Nc electrons, iVn hydrogen atoms, and A^moi 
molecules in a volume V, the number densities being 
Uj = Nj/V. The Helmholtz free energy is written as the 
sum 

= + ^^i'd + FiT + Fg + F^r, (1) 

where i^^, FP^, and are the free energies of ideal 

gases of the electrons, protons, and neutral species, re- 
spectively, F(£ takes into account the Coulomb plasma 
nonideality, and -F^x " nonideal contribution which 

arises from interactions of bound species with each other 
and with the electrons and protons. Ionization equilib- 
rium is given by minimization of F with respect to par- 
ticle numbers under the stoichiometric constraints, keep- 
ing V and the total number A'o of protons (free and 
bound) constant. The latter number is determined by 
the total mass density p: uq = Nq/V ~ p/mn, where 
mn — TTic + TOp. The formulae for each term in Eq. |^ 



are given in IPotekhin et all l)1999() and in Paper I. The 
employed minimiz ation technique i s similar to the tech- 
nique presented bv lPotekhinI l)1996|) for the nonmagnetic 
case. 

Once the free energy is obtained, its derivatives over p 
and T and their combinations provide the other thermo- 
dynamic functions. 

The atomic number fraction xh — n^^/riQ, evalu- 
ated in the course of the free energy minimization, 
can be used in calculations of atmospheric opacities. 
One should take into account that the strong mag- 
netic fiel d affects the po larization properties of radia- 
tion fe. g. . lGinzburellT97(ll) . At photon energies hw much 

higher than (47r/i^e^ne/mo)^''^ « 28.7/9^^^ eV, where 
Pq = p/(l g cm"'^), radiation propagates in the form 
of two so-called normal modes. These modes have dif- 
ferent polarization vectors Sj and different absorption 
and scattering coefficients, which depend on the an- 
gle 9b between the propagation dire ction and B (e.g., 
iKaminker. Pavlov, fc Shibanovlll982|) . The two modes 
interact with each other via scattering. iGnedin fc Pavlovl 
l)1973l) formulated the radiative transfer problem in terms 
of these modes. 

Let the magnetic field be directed along the z-axis. 
At a fixed photon frequency tj, the absorption opacity 
Kj{6b) in each mode j and scattering opacities K^jj,{OB) 
from mode j in to mode j' can be presented as (e.g., 
IKaminker et al.lll982|) 

1 

n^{eB)^m^' J2 |e,.o(0B)P<, (2) 

a = -l 

X— He,,, Sine's de's, (3) 
"^H Jo 

where a — 0, ±1, ej^ — Cj^z is the z-component of e^, 
andcj ii = (ej^a;iiej.j,)/'\/2 are the circular components. 
The cross sections aa depend on u, but not on j or 9b- 
The total scattering opacity from mode j is = K^ji + 
K^2i and the total extinction opacity is Kj = + k^. 
According to equations (0, Q, we can write 

1 

^j{Ob)^ Y.\eUeB)\' k^, (4) 

where Kq {a = —1,0, 1) do not depend on Ob- With a 
good accuracy, ka « (cTq + cr^)/mH. 

In a partially ionized atmosphere, the opacity is con- 
tributed by electrons, ions, and bound species. The scat- 
tering cross section includes contributions from the elec- 
trons and protons: cr^ = + ct^'P. With a good ac- 
curacy, a nd ct^'P are describe d by simple analytical 
formulae fe.g.. lPavlov et anil99lTt) . The absorption cross 
section includes contributions from absorption by 
plasma electrons and protons (mainly by free-free tran- 
sitions due to the electron- proton collisions, a^), transi- 
tions between discrete states of an atom (bound-bound 
absorption, a]^) and photoionization (bound-free, crjj^). 
Thus, for the hydrogen atmosphere, we can write 

<«^H(^^^ + + (l-^H)fTf. (5) 
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For the hydrogen atoms, electrons, and protons mov- 
ing in a strong m agnetic field, the c r oss se ctions a]^ 
were studied by 'Pavlov fc PotekhinI l)1995D . a^' by 
^otekhin fc Pavloy (1997), and cr„ in Paper I. 

Ho & Laii n|200.ll) presented convenient formulae for 
\^j,a{0B)\'^ for a fully ionized, strongly magnetized 
electron-ion plasma (taking into account vacuum polar- 
ization). However, the bound species affect the dielec- 
tric tensor of the medium and hence the polarization 
properties of the normal modes. For a completely neu- 
tral hydrogen gas in a s trong magneti c field, polarization 
modes were studied bv lBulik fc Pavlo v (1996), but for a 
partially ionized gas the problem is not yet solved. For 
the time being, we use the formulae of Ho & Lai ( 200^ . 
Since the neutral fraction is typically small, we expect 
that the resulting error should not be large for the effec- 
tive Rosseland mean opacities presented in the tables. 

3. RESULTS 

3.1. Calculation of tables 

Since our model of the ionization equilibrium and equa- 
tion of state is computationally expensive, it is not pos- 
sible to use our computer code "on line" in any practical 
application. The alternative is to tabulate thermody- 
namic quantities covering the density, temperature, and 
magnetic field domain of interest and to rely on an inter- 
polation of the tabular values. In Paper I, we presented 
the tables which cover a range of p, T, and B appropri- 
ate for most typical neutron stars, such as isolated radio 
pulsars. We used the fitting formulae for the binding 
energies, quantum-mechanical sizes, oscjUator strengths, 
and electron collision widths derived bv lPotekhin (J,99^ . 
Calculations at i? > 10^^-^ G were hampered by the ab- 
sence of the relevant fits. 

In order to overcome this difficulty, in the present work 
we have calculated the tables of the required atomic 
quantities as functions of the transverse pseu domomen- 
tum K I , using the computer code described bv lPotekhinI 
()1994|) . A representative grid of w 50 equally spaced 
values of logK± has been used for every relevant set 
of atomic quantum numbers and every considered field 
strength. Handling of these tables is easier at super- 
strong than at weaker fields, because the number of in- 
volved discrete states is smaller. For example, there 
remain no bound states with s — 2 and s = 1 at 
B > 2 X 10^3 G and B > 6 X 10^^ G, respectively. 
The latter point is illustrated by Fig. ^ which shows 
the energies of the most important atomic transitions at 
K± = as functions of B. With increasing B, the tightly 
bound levels s = 3, 2, and 1 consecutively merge with the 
continuum and cease to contribute to the bound-bound 
spectrum; instead, they appear as autoionization reso- 
nances. There still remain excited loosely bound states 
(with s = and any u) at any B, but, at the typical 
neutron-star atmosphere densities, most of them are de- 
stroyed by interactions with surrounding particles. 

We use a cubic interpolation across these precalculated 
tables in our computer code for obtaining the number 
fractions of atoms and molecules, thermodynamic func- 
tions, and cross sections a]^. As previously, we employ 
the approximate expression for [Eq. (51) in Paper I] 
and numerical tables of crjj^ as functions of K± and oj 
(§4.2 of Paper I). 




9 10 11 12 13 14 15 
log B (G) 

Fig. 1. — Energies of radiative transitions for the H atom with 
K± = from the ground state to the lowest tightly bound states 
{u = 0, s = I, 2, 3; light solid lines) and loosely bound state (s = 0, 
u = I; dashed line); the photoionization threshold (heavy solid 
line); and the energy of the proton cyclotron resonance (dot-dashed 
line). The hatched region corresponds to the bound- free transi- 
tions. 



In this way we have calculated and tabulated ther- 
modynamic functions and opacities of the hydro- 
gen plasma at 13.5 < logS [G] < 15.0 (prelim- 
inary results for B = 10^** G w ere presented by 
iChabrier. Douchin. fc Potekhinll2002|) . We have also up- 
dated the Fortran program for the cubic-polynomial in- 
terpolation of the tabulated data in the 3-parameter 
space of B, T, an d p.^ 

As discussed bv lPotekhin et al.1 1)1999|) . our model be- 
comes less reliable at relatively low T and high p, partic- 
ularly because of formation of molecules and chains H„ , 
which are treated in an approximate manner. In this do- 
main, the partial number fractions and thermodynamic 
quantities are strongly model-dependent. The stronger 
magnetic field, the higher the tempe rature bel ow which 
the molecules and chains dominate ijLail 12001]) . Trying 
to avoid the domain of uncertainty, we shifted the lower 
bound on the temperature in the tables, as compared to 
the case of weaker fields (Paper I), from 2 x 10^ K to 
5 X 10^ K. On the other hand, the temperature at which 
the atoms become thermally ionized increases with in- 
creasing B, therefore we have shifted the upper tempera- 
ture bound in the tables from 10'^ K to 4x 10^ K. It would 
be meaningless to consider still higher T, because hydro- 
gen undergoes efficient the rmonuclear burning at such 
temperatures l)Ergmalll986j) . The third input parameter 
in our tables, in addition to B and T, is the "astrophysi- 
cal density parameter" R = po/Tg , where Tq = T/10^ K. 
In order to avoid the rather uncertain region of possible 
phase separation (see below) , we restricted the public ta- 
bles by logi? < 3.4, although higher R values were also 
considered. This range of R should be sufficient for mod- 
eling the atmospheres of most neutron stars, except the 
unusually cold ones mentioned below. 

Along with the pressure, our tables contain internal en- 

^ T he tables and the program are available at 
http://www. iof f e . ru/astro/NSG/Hmagnet/ 
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Fig. 2. — Neutral fraction of the H atoms in all quantum states (heavy lines) and in the ground state (dashed lines) as function of density 
at the magnetic field strength B = 0, 10^^'^ G, lO^"^'^ G, and 10^^'^ G. The dot-dashed line shows the fraction of H2 molecules at the 
highest B. Left panel: T = 5 X 10^ K; right panel: T = 2 X 10*^ K. 
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Fig. 3. — Domains of partial ionization at logB [G] = 14.0, 14.5, and 15.0. The contours delimit the domains where the atomic fraction 
exceeds 0.1% (dotted lines), 1% (dashed lines), or iO% (solid lines). Hatched is the domain where the H2 fraction exceeds 1%. 



ergy, entropy, specific heat, and the logarithmic deriva- 
tives of pressure with respect to V and T. Other second- 
order thermodynamic quantities can be ca lculated using 
the Maxwell relations ijLandau fc Lifshi't3[l993y) . 

The format of the tables is the same as in Paper I. 
In addition to the thermodynamic quantities, the ta- 
bles contain the number fractions of different chemical 
species, and also the effective Rosseland mean opacities 
for diffusion of radiation along and across the magnetic 
field. 

3.2. Ionization equilibrium and equation of state 

A strong magnetic field generally increases the frac- 
tion of bound species at a given T. Figure |2 illustrates 
this point. In this figure, ionization equilibrium curves 
are plotted as functions of density for different values 
of B and T. At the lowest T and highest B shown in 
Fig. 121 the molecular II2 fraction (dot-dashed line on the 
left panel) becomes significant at p > 1 g cm^^. At the 
same T and B and higher density p ~ 100 g cm^"^, the 
II2 molecules disappear, but the II„ chains and clusters 



(not shown in the figure) become abundant. As men- 
tioned above, our results are not well justified under such 
conditions. 

Figure 13 shows the domains of partial ionization in the 
T-p plane at three values of B. With increasing B, the 
domains where the atomic fraction is above a specified 
level expand significantly towards higher temperatures. 

Figure0]shows pressure P as function of p at T = 10®'^ 
K and at different field strengths. At very low p, we have 
nearly fully ionized, almost ideal nondegenerate gas of 
electrons and protons. The magnetic field does not af- 
fect the equation of state in this case. With increasing 
density, at lO'^ g cm^'^ ^ P ^ 10'^ g cm~'^, the electrons 
become degenerate at i? = 0, but remain nondegenerate 
at B > 10^^ G. This explains why the dashed line in the 
figure (pressure at i? = 0) is higher than the other lines 
in this regime. Meanwhile, because of the partial recom- 
bination of atoms, the solid curves {P of the nonideal 
magnetized plasma) are slightly lower than the dotted 
ones (P of the ideal fully ionized plasma) . At still higher 
density, however, the atoms become pressure ionized: in 
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Fig. 4. — Pressure isotherms (T = lO'''^ K) of partially ionized 
hydrogen in the strong magnetic field (B = 10^'^ G, 10^* G, and 
10^^ G; solid lines) compared with the model of a fully ionized 
ideal electron-proton plasma (dotted lines) and with the case of a 
hydrogen plasma at B = (dashed line). The vertical lines cor- 
respond to the density above which excited Landau levels become 
populated. The instability region at _B = 10^^ G (the discontinu- 
ity of the lower solid curve) corresponds to a hypothetical plasma 
phase transition (see text). 



this region the increase of P due to the increased num- 
ber of free electrons and protons competes with a neg- 
ative nonideal contribution, which is mainly due to the 
Coulomb term in the free energy. 

At sufficiently high B or low T, the Coulomb interac- 
tion leads to a discontinuity, exemplified in Fig. 01 by 
the line B = 10^^ G. It occurs at T below the crit- 
ical temperature Tc and a t p around the critic al den- 
sity pc ~ 143^1^2^8 g cni-3 l)Potekhin et a l."1999'). where 
Bi2 = B/10^^ G. This phase transition is directly related 
to the magnetic field assisted separation b etween the con- 
densed hydrogen and vapor, discussed bv lLai fc Salpete^ 
(^997). Note that the density of the condensed phase 
at zero pressure and zero temperature, estimated by 

iLai & Salneie^ l|1997f) as ps « 560 5^^^ g cm'^, scales 
with B similarly to pc- At B = 10^^ G, this estimate 
gives Ps ~ 2.2 X 10^ g cm~'^, whereas we have numeri- 
cally Ps ~ 1-7 X 10^ g cm~^ in Fig. 21 However, in gen- 
eral, the validity of the free-energy models in the frame- 
work of the chemical picture of plasmas is questionable 
near the plasma phase transition. In order to study the 
matter properties near the phase transition, one needs 
to consider accurately the interaction of H„ molecules 
in the gas phase and their interaction with the surface 
of the condensed phase. This problem goes beyond the 
scope of the present paper, although it certainly deserves 
a future study. If this phase separation is real, it can 
be important for interpretation of the emission of neu- 
tron stars with unusually low T and/or high B, which 
in this case do not possess an optically thick atmosphere 
above the condensed surf a ce. T his hypothesis was sug- 
gested by ILai fc Salpeteil l)1997j) ; a possible example of 
su ch "naked neu t ron st ar" has been recently considered 
bv lTurolla eTal] l|2nf)3|) . 



3.3. Opacities 

Figure |31 shows the monochromatic extinction coeffi- 
cients ka = pka for the three basic polarizations a — 
0,±1 aX B = 10^3.5 G, Tg = 1 or 4, and p = 0.1 g cm-^ 
or 50 g cm^'^. Figure demonstrates the same func- 
tions for 10 times higher magnetic field and densities. 
These densities correspond to the outer and inner lay- 
ers of the atmosphere, where the optical depth is of the 
order of 1 for ordinary (high-opacity) and extraordinary 
(low-opacity) rn odes of radiation, respectively (see, e.g., 
iHo fc Laill2003t) . At the higher densities, the opacities 
are featureless, except the proton cyclotron resonance at 
hijJc-p = hijJcciTic/m-p — 6.305 i?i2 eV, because the plasma 
at such density is pressure-ionized. On the contrary, ab- 
sorption features are clearly visible at the lower density. 

At the lower field strength (Fig.EJl, one sees the broad 
absorption feature for a = ±1 (the left and central pan- 
els) due to the s 1 bound-bound transition at 
hio < 0.3 keV. For a = -1-1 (the left panel), this ab- 
sorption merges with the strong proton cyclotron free- 
free resonance at lower energies and with the photoion- 
ization continuum at higher energies. This situation 
is rather special, because at this field strength hiOcp is 
rather close to the photoionization threshold, whereas 
the limiting energy for the bound-bound transition falls 
in between (see Fig. . Both the photoionization edge 
and the absorption line are considerably broadened be- 
cause of the energy shifts du e to the a tomic motion 
(the magnetic broadening: se e iPavlov fc~P otckhin 192^ 
and lPotekhin fc Pavlovlll997j) . therefore they effectively 
merge together. On the other hand, the free-free absorp- 
tion coefficient at the lower T (the solid curve) reveals 
the bumps at ^iw < 1 keV, which arise from the opening 
of ionization channels with higher s values at higher ener- 
gies, and from autoionization resonances just below thes e 
partial ionization thresholds l)Potekhin fc Pavlovlll997f) . 
On the right panel (a = 0, the polarization vector along 
B), one notices only the photoionization edge at rela- 
tively small p and T (the lower solid line). 

At the higher field strength (Fig. EJ, all transitions 
changing s belong to the continuum (see Fig.^ and form 
a series of resonances on the photoionization profile at 
the energies above fi^cp, particularly visible for a = -1-1 
at low p and T (the lower solid line on the left panel). 
Near the bottom of the photosphere (p = 500 g cm"'^), 
the bound states are destroyed and do not contribute to 
the spectrum. 

4. CONCLUSIONS AND OUTLOOK 

We have calculated the equation of state of fully and 
partially ionized hydrogen plasmas in a wide range of 
densities, temperatures, and magnetic fields expected for 
photospheres of the magnetars. We have also calculated 
monochromatic radiative opacities for three basic po- 
larizations. The first- and second-order thermodynamic 
functions, non-ionized fractions, and effective Rosseland 
mean opacities are published in the electronic form. 

These results can be useful for construction of the at- 
mosphere models and calculation of the spectra of out- 
going radiation from warm neutron stars with super- 
strong magnetic fields. This work is under way ; some 
results have been published bv lHo et al.l l)2003all]|) . How- 
ever, there remain unsolved problems which hamper an 
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Fig. 5. — Monochromatic extinction coefficients for the three basic polarizations, = p^ai = 0, ±1, at p = 0.1 g cm or 50 g cm ^, 
T = lO" K (sohd hnes) or 4 x lO" K (dashed hnes), B = lO^^-S q. 
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The same as in Fig.|21 but for p = 1 g cm~'^ or 500 g cm~'^, B = 10^*'^ G. 
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immediate application of our present results to model- 
ing accurate and reliable outgoing spectra of magnetars 
with hydrogen atmospheres. First, vacuum polarization 
alters the dielectric tensor of the medium and polar- 
ization of photon modes (see iPavlov &: GnedinI [l98^ . 
which is important for magnetar atmospheres. In par- 
ticular, the "vacuum resonance" phenomenon may lead, 
under certain conditions, to an adiabatic mode conver- 



sion and "mode collapse" ijHo fc Lal2003l and references 
therein). In the latter case, the description in terms 
of the two modes fails, and one has to solve more gen- 
eral equati ons for the evolution of electromagnetic waves 
()Lai fc Hdl2003|) . Second, as noted in 32 the presence of 
the bound species also changes the dielectric properties 
of the medium and the normal mode polarization, which 
necessitates derivation of a polarization tensor for par- 
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tially ionized, strongly magnetized plasmas. Third, for 
relatively cold neutron stars with strong magnetic fields, 
the phase separation mentioned in t |3.2l should be studied 
in a thermodynamically consistent way. We are planning 
to study these problems in near future. 
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